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Abstract 

Recent measurements of time dependent CP asymmetry in B — > cpKs, if confirmed, 
would indicate a new source of CP violation. We examine flavor violating tree-level Z 
currents in models with extra down-type quark singlets that arise naturally in string 
compactified gauge groups like Eq. We evaluate the new operators at the scale /i « 
0{rrib) in NLO, and using QCD improved factorization to describe B — > 4>Ks, find 
the allowed range of parameters for p and tjj, the magnitude and phase of the flavor 
violating parameter Ub s - This allowed range does satisfy the constraint from flavor 
changing process b — ► s£ + £~ . However, further improvement in measurement of these 
rates could severely constrain the model. 
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1 Introduction 



The ongoing B physics experiments by BaBar and Belle collaborations [1, 2] provide a 
unique opportunity to study the flavor structure of the standard model quark sector and 
also the origin of CP violation. In addition to this, any new physics effects in B physics 
can also be tested in these experiments. Recent time dependent asymmetries measured in 
the decay B — > (j)K s both by BaBar and Belle collaborations [1-4] show significant deviation 
from the standard model and this has generated much theoretical speculation regarding 
physics beyond the standard model [5]. In the standard model, the process B — > 4>K S is 
purely penguin dominated and the leading contribution has no weak phase. The coefficient 
of sin(Amst) in the asymmetry therefore should measure sin 2/3, the same quantity that is 
involved in B — > ipK s in the standard model. The most recent measured average values 
of asymmetries are [4, 6] S^k s = 0.734 ± 0.055 and S<pK s = —0.15 ± 0.33. The value for 
SipK s agrees with the standard model expectation. The deviation in the <t>K s is intriguing 
because a penguin process being a loop induced process is particularly sensitive to new 
physics which can manifest itself through exchange of heavy particles. In this article we will 
consider an extension of the standard model, with extra down type singlet quarks. These 
extra down type singlet quarks appear naturally in each 27-plet fermion generation of Eq 
Grand Unification Theories (GUTs) [7-10]. The mixing of these singlet quarks with with the 
three SM down type quarks, provides a framework to study the deviations from the unitarity 
constraints of 3 x 3 CKM matrix. This model has been previously studied in connection 
with Rb and F-B asymmetry at the Z pole as it provides a framework for violation of the 
unitarity of the CKM matrix [10-12]. This mixing also induces tree-level flavor changing 
neutral currents (FCNC). These tree-level FCNC couplings can have a significant effects on 
different CP conserving as well as CP violating B processes [11, 13 - 24]. 

In this article we study the FCNC effect arising from the Z — b — s coupling Ub s to the 
B — > 4>K S process. This new FCNC coupling U bs can have a phase, which can generate the 
additional source of CP violation in the B — > 4>K S process, and thus affect measured values 
of S(f,K s and C^,k s - We parameterize this coupling by \J hs = pe 1 ^. We then study B — > 4>K S 
taking into account the new interactions in the QCD improved factorization scheme (BBNS 
approach ) [25]. This method incorporates elements of naive factorization approach (as its 
leading term ) and perturbative QCD corrections (as sub-leading contributions) and allows 
one to compute systematic radiative corrections to the naive factorization for the hadronic 
B decays. Recently, several studies of B — > PV, and specifically B — > 4>K S have been 
performed within the frame work of QCD improved factorization scheme [26-30]. In our 
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analysis of B — > (f>Ks, we follow [30] which is based on the original paper [25]. In our analysis, 
we only consider the contribution of the leading twist meson wave functions, and also neglect 
the weak annihilation contribution which is expected to be small. Inclusion of these would 
introduce more model dependence in the calculation through the parameterization of an 
integral, which is otherwise infrared divergent. 

The time dependent CP asymmetry of B — > 4>K S is described by : 

A m rgog) -> <j>K s ) -T(B°(t) -> <j>K s ) 

* Ks[) ~ r (W(t) - 4>k s ) + r(B°(t) - 4>k s ) [i) 

= -C^Ks cos(Am B t) + S^Ks sin(Am B t) (2) 
where S^Ks an d ^Ws are given by 
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b <t>K S - , i x [a > - 7T I \ 12 ^ 

J-+ I *4>k s | J-+ | -Vi^s I 

and X(/>k s can be expressed in terms of decay amplitudes: 

_ M(B^<PK S ) 

X * Ks ~ " 6 M{B^Ks) (4) 

The branching ratio and the direct CP asymmetries of both the charged and neutral 
modes of B — > have been measured [1, 2, 3, 4, 6, 31] 3 

B(B° -> 0AT s ) = (8.0 ± 1.3) x 10~ 6 (5) 

B(B + ^(f)K + ) = (9.4 ±0.9) x 10~ 6 , (6) 

S 4>Ks = +0.45 ±0.43 ±0.07 (BaBar); (7) 

= -0.96±0.50±°;?i (Belle); (8) 

= -0.15 ±0.33 (World average); (9) 

C 4>Ks = -0.19 ±0.30 (10) 

A C p(B + -> (f)K + ) = (3.9 ±8.8 ±1.1)% (11) 



3 Latest results were reported at XXXIX Rencontres de Moriond, Electroweak Interactions and Unified 
Theories, Italy, March 2004. See talks [32, 33, 34]. The Belle result on S^k is unchanged [32] while BaBar 
finds Sqk — 0.47±0.34±°;°8 [33] which is very close to the result in Eq.(7). Hence, our observations remain 
unchanged. 
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scale 


Ci 


c 2 


c 3 


C 4 


c 5 


c 6 


fj, = m b /2 
H = m b 


1.137 
1.081 


-0.295 
-0.190 


0.021 
0.014 


-0.051 
-0.036 


0.010 
0.009 


-0.065 
-0.042 




CV/ «era 


C$/ a em 


C%j a em 


Clo/ «em 




C 8 g 


fi = m b /2 


-0.024 


0.096 


-1.325 


0.331 


-0.364 


-0.169 


H = m b 


-0.011 


0.060 


-1.254 


0.223 


-0.318 


-0.151 



Table 1: Standard model Wilson coefficients in NDR scheme. 



2 B — > in the QCDF Approach 

In the standard model, the effective Hamiltonian for charmless B — > 0i^5 decay is given by 
[25] 



^eff = ~^ VtbV ts 



cmom + cmom + 1 £c i {ti)o i (ti) + c 77 o 77 + c 8g o 8g 



10 

(12) 

i=3 

where the Wilson coefficients Cj(//) are obtained from the weak scale down to scale /j, by 
running the renormalization group equations. The definitions of the operators can be found 
in Ref.[25]. The Wilson coefficients Cj can be computed using different schemes [35]. In this 
paper we will use the NDR scheme. The NLO values of Cj(i = 1 — 10) and LO values of 
£77; C 8g respectively at /i — m b /2 and m b used by us based on Ref. [25] are shown in Table 1. 

In the QCD improved factorization scheme, the B — > (pKs decay amplitude due to a 
particular operator can be represented in following form : 

<(j)K\0\B >=< <t>K\<D\B> fact [l + Yl r n< + 0(A QCD /m b )] (13) 

where < 4>K \ \ B >/ act denotes the naive factorization result. The second and third term 
in the bracket represent higher order a s and AQ C D/ m b correction to the hadronic transition 
amplitude. Following the scheme and notations presented in Ref. [30], we write down the 
total B — > (f)Ks amplitude, which is the sum of the standard model as well as Z exchange 
tree-level contribution from extra down-type quark singlets (EDQS) model in the heavy 
quark limit 

M{B + ^<PK + ) = 

M(B° — > 4>K°) = ^mlUFr K H)V pb V; s [al + < + < 



^2 ""&J<P~ 1 V"(p/ • PO • ps 

P + nP + n p „\ 1 

(14) 



(a? + ag + a? ) . P 

o + «10a 
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where p is summed over u and c. The coefficients af are given by 



C' 4 

Go = do = Cq H — 



4:71 



a' 

a 4 - G 4 + 



4tt 



C F a s 
4jrN r 



PI 



< = at = C' 5 + f 



CL'j — Cly — C<j ~h 



<5 

N 

a 



c 

I r 

8 



n u _ n c _ r i , °10 
a 9 — a 9 — °9 + 



1 + 
1 + 



47T 
CfCXs 
47T 

Cpa t 



ci-\ rt — a 



V 

aio 



10 ~~ "10 ~~ 



Cl ° + TV, 



1 + 



47T 



47T 



-12 -V^-H^) 
(V* + if,) 



4tt^ 



Ql 



a 10a — a 10a 



CpOL t 

4txN ( 



Q<t> 



(15) 



with C F = {N 2 - 1) /2N C and iV c = 3. The effective Wilson coefficients C[ = Q+Q, (i = 3- 
10) is the sum of the standard model and the EDQS model Wilson coefficients. The quantities 
Vfj), H^, and are hadronic parameters that contain all nonperturbative dynamics, are 
given in Ref. [25, 36]. 

For the sake of completeness, we give the branching ratio for B — > <f)K,s decay channel in 
the rest frame of the B meson. 



BTZ{B -> (j)K s ) = 



TB I Pc , 



M(B^<PK S ) 



8n m 2 B 

where, r B represents the B meson lifetime and the kinematical factor | P cm | is written as 

, _ , 1 



(16) 



2m i 



m 2 B - (m K + m^) 2 } [m 2 B - (m K - m^) 2 ] 



(17) 



3 FCNC Z couplings in EDQS model 

Models with extra down-type quarks (EDQS) have a long history. The earliest consideration 
of such models was in the context of the grand unification group E$ which arises from string 
compactification. The quarks and leptons of each generation belong to the 27 representation 
[7-10]. Each generation has one extra quark singlet of the down type, and also one extra 
lepton of the electron type. The group also has extra Z bosons, which we will assume to 
be too heavy to have any effect on B — > <pKs process. The down type mass matrix is then 
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a 6 x 6 by-unitary matrix, and in general when we rotate the quarks to their mass basis, 
off-diagonal couplings arise. In EDQS model, the Z mediated FCNC interactions are given 
by 

c = ^hr [dLaU aP rd L p] z, (is) 

z cos t>w 

In general for n copies of extra down-type quark singlet model, U a p is : 

3 N d 

U a p = E VliVifi = <W - E v LViP, («> P = d,s, b, B u B 2 ....) (19) 

i=l i=4 

where, N d — 3 + n represents the number of down type quark states, and U is the neutral 
current mixing matrix for the down quark sector. The non vanishing components of U a p will 
lead to FCNC process at tree level, generating new physics contribution to the measured CP 
asymmetries. The new tree level FCNC Z mediated contribution to the b — > sqq process is 
shown in Figure 1. The new operators arising from this tree- level FCNC process have been 
shown to lead to the following effective Hamiltonian for b — > sqq process in this model [22]: 



ri 7 = 7= Vth V+ 



G F 



C 3 3 + C 7 7 + C 9 9 (20) 



where, the new Wilson coefficients C 3 , C 7 and C 9 at the scale M z are given by: 

C 3 (M Z ) = (21) 

C 7 (M Z ) = K ^sm 2 9 w , (22) 

C 9 (M Z ) = -^(l-sin 2 ^) (23) 

where, k = ) , and operators Oi in Eq. (20) are given in Ref.[25]. We now evolve these 
new Wilson coefficients from the scale M z to the scale \x 0(m b ) using the renormalization 
group equation. While doing this we have considered NLO QCD correction [37], neglecting 
the order a electroweak contributions to the RG evolution equation which are tiny. At the 
low energy, after the RG evolution the above three Wilson coefficients (C3, C 7 , C 9 ) generate 
new set of Wilson coefficients (Ci,i = 3 — 10) in this model. The values of Wilson coefficients 
( without taking the overall factor k) at scales ji = (mb/2,mb) are shown in Table 2. 



4 B physics constraints on Ub s 

In this section, we review the constraints on the flavor violating parameter Ub s from different 
flavor changing B processes. These processes can be classified into two classes, CP conserving 
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scale 


c 3 


c 4 


c 5 


c 6 


C 7 




c 9 


Cio 


fj, = m b /2 


0.195 


-0.088 


0.0180 


-0.053 


0.133 


0.108 


-0.604 


0.174 


fj, = mb 


0.182 


-0.0629 


0.0157 


-0.0370 


0.136 


0.0732 


-0.574 


0.122 



Table 2: Wilson coefficients of EDQS model in NDR scheme, without the overall multiplica- 
tive factor k. 



and CP violating. Among the CP conserving processes, B(B — > X s £ + £ ), and AM Ba can put 
constraints on U bs [13, 15, 16]. Using recent Belle [38] measurement of B(B — > X s l + l~) = 
(6.1 ± 1.4±i;i) x 10" 6 the authors in Ref.[21] had shown that | U bs |< 1 x 10~ 3 . However, 
this bound has recently been updated in Ref. [22] to 

| U bs - 4.0 x 10~ 4 |< 8 x 10~ 4 (24) 

which also updates their previous bounds in Ref. [17, 39] The bound in Eq. (24) is based on 
inclusive B — > X s e + e~ decays at NNLO [40]. We shall adopt this bound in our analysis. 
This bound is valid in both general n extra down-type singlet quark model and in the model 
with a single extra down-type quark singlet [16]. Similarly, the b —> s'j branching ratio 
provides comparable limits [14, 16]. 

It has been shown in Ref. [20, 41] that in the presence of tree-level FCNC coupling Um 
and /or U bs , the standard box diagram for B d / S — B d / S is not gauge invariant by itself, but 
requires Z exchange penguin diagram as well as tree level Z FCNC diagram. The additional 
Feynman diagrams are given in Ref. [41]. Following the paper [20], with slight change in 
the notations, we write down the expression for the B q — B q mixing: 



AMn = 



G 2 F M^B Bq m B0q 



6tt 2 

where, (q = d, s) and 



(Xl b )\f t s S (x t ) + A r , 



47rsin 2 #w B 2 



(25) 



A new = -8U bq \l b r)? t <>Y (x t ) + -Vz q Ui (26) 



a 



where, the definitions of different parameters used above can be found in Ref. [ 

We have found that to satisfy the measured AM Bd 4 within one sigma, where we consider 
both the theory error of 20% arising from the value of f% B Bq and the experimental error 



*The experimental value for | AM Bd \= 0.489 ± 0.008 ps" 1 [20]. 
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taken in quadrature, the FCNC coupling | U bd | should be less than ~ (2 — 3) x 1CT 4 . This 
is a very stringent limit. 

The AM Bs has not been measured yet, and so only lower limit on the mass difference 
is available. We have found that the new contribution to AM Bs from EDQS model is less 
than 3% when compared to the standard model. It can be shown that for similar values of 
U bd and U bs , the FCNC effects on AM Bs will be suppressed by a factor ~ A 2 when compared 
with the effects on AM Bd . This implies that in EDQS model, FCNC effects are hard to 
detect in B s — B s mixing [41]. 

5 B — » cj)Ks analysis 

In the last section we have discussed the allowed range of the FCNC parameter U bs from 
different B processes. In this section we will study the effect of this FCNC parameter (U bs ) 
in the B — > 4>Ks process. For this we express U bs in the following form: U bs = pe 1 ^. We will 
then vary p and ip 5 in range such that Eq. (24) is satisfied. We then study the allowed region 
of parameters in the p — ip plane from the three measured quantities (a) E[B — > (f)Ks), (b) 
S<pK s and (c) C^k s - To get the allowed parameter space, from the B(B — > <pK s ) branching 
ratio, we allow it to vary by 2a respectively from its central value. This 2a band contains 
both experimental and theoretical errors. The main source of theoretical error is the form 
factor Ff^ K . In our analysis we have considered 20% error on this parameter. Similarly, 
we vary C^k s and S^Ks by la and 2a from their central value to get the allowed region in 
the p — ip plane. 

In Figure 2 (a) we show such allowed region in p — ip plane for the scale p = m b /2. 
The whole area left of the dotted contour is allowed by saturating Eq. (24). The area 
outside the thick contour labeled by BR is 2a allowed region from the branching ratio 
measurement. The parameter space enclosed by the thin contour marked by S^k is allowed 
by 2a from data on the S^k s - This whole parameter space is allowed by la from C^k s 
measurement, The regions (marked by Z) is the only allowed parameter space in p — ip plane 
with 6.5 x 10~ 4 < p < 10 x 10~ 4 and —1.7 < ip < —0.85 which satisfy the experimentally 
measured C^Ksi S<t>K s and B(B — > <pKs) within the errors described above. We note that 
only negative values of ip give acceptable range of S^k- The Figure 2 (b) correspond to the 
scale p = m b . In this case though we have larger allowed area from the S^k measurement, 
but the 2a branching ratio contour pushes the allowed range towards higher values of p and 

5 ?/> in units of radian 
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somewhat lower range of the phase ip. This particular behavior of the branching ratio contour 
can be understood from that fact that for p = m b , the SM branching ratio is 3.8 x 1CT 6 , 
which is much smaller than the lower end of the 2a band of the experimental number. Hence, 
one needs larger values of p to push the total branching ratio within the 2a limit. For this 
reason the allowed region shrinks to a point in this case. 

6 Conclusions 

In this paper, we have studied the tree-level flavor violating Z contribution to B — > <pKs 
process in models with extra down-type quark singlets which arise naturally in the context 
of the Grand Unification group Eq. In the presence of such flavor violating interactions, 
B — > <fiKs process receives additional contributions, governed by a set of new operators which 
can be expressed in terms of the standard operators Oj, (i = 3 — 10). We then evolved these 
new Wilson coefficients from the scale Mz to the scale p = 0{m^) relevant for our process 
using the renormalization group equation. We have found that, at the lower scale these 
Wilson coefficients significantly modified from their initial values at the scale p = M z . We 
have found that this new flavor violating interaction can modify the standard model Wilson 
coefficients C i: (i — 3 — 10) significantly. We have used following experimentally measured 
quantities: S^Ks^C^Ks an d B(B — > <pK s ) to constrain the flavor violating parameter [/& s . 
We have shown that in the model with an arbitrary number of down type singlet quarks, 
the value of S<j,K s an d C^k s can be well explained by the values of p and tp in the region 
marked by Z in Figure 2 (a). Improvements in measurements of B — > X s £ + £~ can tighten 
the constraints in Eq. (24) and either rule in or rule out this model. 
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7 APPENDIX : Input parameters and different form 
factors 

In this Appendix we list all the input parameters, decay constants and form factors used for 
the calculation of B — > 4>K S . 

1. Coupling constants and masses ( in units of GeV ): 

a em = 1/129, a s (M z ) = 0.118, Gf = 1.16639 x 1CT 5 (GeV) -2 , 
M z = 91.19, m b = 4.88, m B = 5.2787, 
= 1.019, m K = 0.493 

2. Wolfenstein parameters : 

A = 0.2205, A = 0.815, ?? = 0.324, p = 0.224, 

3. Constituent quark masses rriiii = u, d, s, c, d)( in units of GeV): 
m u = 0.2, m d = 0.2, m s = 0.5, m c = 1.5, m h = 4.88. 

4. The decay constants (in units of GeV): 

/b = o.i9, iV = 0.237, f K = o.ie 

5. The form factors at zero momentum transfer : 
F^ K = 0.33. 
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Figure 1: Feynman diagram for Z exchange tree-level contribution to b — > sss process. 
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Figure 2: Contour plots of S^Ks, B(B — > <pK s ) and | U bs -A x 1(T 4 |= 8 x 10~ 4 in p — tj; plane 
for two values of /i — ^ (a) and m b (6) respectively. The area leftside of the dotted contour 
is the allowed region of U bs from the inclusive B — > X s £ + £~ process. In figure (a), the area 
marked by Z is the 2a allowed regions from the measurement of S^k and B(B — > 4>K S ). In 
figure (6), such 2a allowed region is a point where the three curves intersect. 
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